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We give calculated results of differential attenuation and differential 
phase shift caused by rain, based on scattering of a plane electromag- 
netic wave by horizontally disposed oblate spheroidal raindrops. Two 
polarizations of the incident wave are considered, as depicted in Fig. 1. 
The factor e~ iut has been suppressed. Subscripts I and II designate 
electric fields parallel and perpendicular to the plane containing the 
axis of symmetry of the raindrop and the direction of propagation of 
the incident wave, respectively. The angle between the direction of 
propagation and the axis of symmetry is denoted by a. For terrestrial 
microwave relay systems we take a = 7r/2, but for satellite systems 
other values of a are of interest. 

The incident wave induces a transmitted field in the interior of the 
raindrop, and a scattered field. In the far field the quantities of pri- 
mary interest are the complex forward scattering functions 1 Si (0) and 
<Sn(0). For the two polarizations considered, the polarization of the 
far scattered field is the same as that of the incident wave. However, 
when oblate raindrops are canted with respect to the vertical, the ver- 
tical and horizontal polarizations commonly used in radio communica- 
tions systems will be neither of the two aforementioned polarizations. 
Then cross-polarization coupling will occur as a result of the difference 
between Si(0) and Sn(0). 

It is assumed that the ratio of minor to major axis of the oblate 
spheroidal raindrop depends linearly on the radius a (in cm) of the 
equivolumic spherical drop; specifically a/b = (1 — a). This relation- 
ship is similar to that used by Oguchi. 2 The Laws and Parsons drop- 
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Fig. 1 — Two polarizations of the incident wave. 



size distribution and Best's raindrop terminal velocity are used. 3 Thus, 
for rain rates up to 150 mm/hour, there are 14 different drop sizes, 
a = 0.025 (0.025) 0.35, to be considered. At 20°C, the refractive indices 
# = 5.581 + 2.848i at 30 GHz and N = 6.859 + 2.716i at 18.1 GHz 
were obtained from an elaborate fitting equation in a recently published 
survey 4 of available measured data. Since the calculations at 4 GHz 
were made at an earlier date, the value N = 8.77 + 0.915% was taken 
from the older literature. 

Oguchi 2 has developed a perturbation theory for scattering from 
slightly eccentric spheroidal raindrops, but the first-order approxima- 
tion is expected to be inaccurate for the larger raindrops, which have 
larger eccentricity. Consequently, we have used a matching theory to 
obtain approximate nonperturbative solutions to the problem. Al- 
though the calculations reported on here are for oblate spheroidal rain- 
drops, the procedure may be used for axisymmetric raindrops which 
are not too nonspherical. Full details of the analytical and numerical 
procedures, together with tables of the computed values of Si(0) and 
<Sn(0) for each drop size, will be given at a later date, but we outline 
the approach here. 

Spherical coordinates (r, 6, <p) are chosen with polar axis along 
the axis of symmetry of the raindrop, and origin at the center, as 
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in Fig. 1. The scattered electric field is expanded in the form 
2S(a m „M m „ + bmnN m „), in terms of solutions of the vector wave equa- 
tion 5 satisfying the radiation condition. An analogous expansion is 
assumed for the transmitted field, in terms of vector wave functions 
which are finite at the origin. The complex coefficients in the expansions 
are determined approximately by satisfying the boundary conditions, 
namely the continuity of the tangential components of the total electric 
and magnetic fields across the surface of the raindrop, by a fitting pro- 
cedure. The index in denotes the order of the harmonic in the azimuthal 
angle <p. The incident field may be expressed as a Fourier series in <p, 
with coefficients depending on r and 6. Because of this, and the axial 
symmetry of the raindrop, the problem can be decomposed and the 
boundary fitting carried out independently for each m. 

A large number of points is chosen on the curve which is the inter- 
section of the boundary of the raindrop with the half-plane <p = 0, 
^ 6 ^ 7r, and the coefficients are determined by requiring the bound- 
ary conditions to be satisfied at these points in the least squares sense. 
The advantage of using least squares fitting rather than collocation 
(in which the number of fitting points is equal to the number of un- 
known coefficients, which are then determined by solving a system of 
simultaneous linear equations) is that the boundary conditions are 
satisfied more accurately. At least twice as many fitting points as un- 
known coefficients were used in our calculations. The number of terms 
required to adequately satisfy the boundary conditions depends on 
both the frequency and drop size. For example, at 30 GHz it was neces- 
sary to take max m = 8 and max n = 23. In order to ensure the ac- 
curacy of the leading terms, it is necessary to take more terms in n 
than are really needed in the calculation of the far-field quantities such 
as £i(0) and Su(0). A convergence test was carried out for each drop 
size, by increasing the upper limit of n in the sums by 2 and by 4. In 
a number of cases more than half the capacity of a Honeywell 6070 
computer was used. 

The rain-induced attenuation and phase shift are obtained from the 
forward scattering functions as follows '} 

Aim = 0.434 — £ ReSi.n(0)3l(fi) dB /km 

7T 

$i.ii = -36 — £ImSi.n(0)9l(a)*Vkm 

4?r 2 
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Fig. 2 — Rain-induced differential attenuation. 



where X is the wavelength in centimeters, 91(a) is the number of rain- 
drops with mean radius a per cubic meter, and the summation is taken 
over all drop sizes. For a = 90 degrees, the attenuation and phase shift 
have been calculated for various rain rates at 4, 18.1, and 30 GHz for 
polarizations I and II. The cases a = 50 degrees and 70 degrees have 
been calculated at 30 GHz. The differential attenuation An — Ai and 
the differential phase shift $n — $i are summarized in Figs. 2 and 3. 
Full numerical results will be presented later. 

At 4 GHz the differential attenuation is negligibly small (0.036 
dB/km at 100 mm/hr), and therefore is not plotted in Fig. 2. How- 
ever, there is a differential phase shift of 5 deg/km at 100 mm/hr which 
can induce cross polarization as large as — 10 dB over a long path of 
heavy rain ; indeed, such a value was observed experimentally on one 
occasion by W. T. Barnett. 6 

The calculated differential attenuations for a = 90 degrees in Fig. 2 
are considerably larger than those obtained by first-order perturbation 
theory. Taking into account the canting angle distribution, 7 these 
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Fig. 3 — Rain-induced differential phase shift. 
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larger values of differential attenuation are needed to explain the mea- 
sured difference in attenuations for horizontal and vertical polariza- 
tions. 6 ' 8 

At 18.1 GHz the differential attenuation and the differential phase 
shift will be comparable contributors to rain-induced depolarization. 
However, the differential phase shift remains small for heavy rain rates 
at 30 GHz, and hence the differential attenuation will be the dominat- 
ing cause of rain-induced depolarization at this frequency and pre- 
sumably also at higher frequencies. When a decreases from 90 degrees, 
both the differential attenuation and the differential phase shift be- 
come smaller as intuitively expected. 

The authors are indebted to D. C. Hogg for bringing this problem to 
their attention, to J. McKenna and N. L. Schryer for suggesting the 
matching and least squares fitting approaches, and for several helpful 
discussions in relation to these, to P. A. Businger whose least squares 
fitting subroutine was incorporated into the main program, to Mary 
Ann Gatto who took over the burdensome task of running the main 
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program, to Susan Hoffberg who wrote the program for calculating the 
first-order approximation for small eccentricities, and to Diane Vitello 
who performed the summation over the drop size distribution. The 
authors are particularly indebted to J. McKenna and D. C. Hogg for 
their continued encouragement throughout the lengthy course of this 
work. 

Note added in proof: 

A very recent paper by Oguchi, 9 in which similar calculations are 
carried out for a = w/2 at 19.3 and 34.8 GHz, has come to our at- 
tention. He used collocation for the expansions in terms of spherical 
vector wave functions, and at 34.8 GHz he also used an expansion in 
terms of spheroidal wave functions. 
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